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ABSTRACT 

Ram pressure stripping can remove significant amounts of gas from galaxies in clusters and massive 
groups, and thus has a large impact on the evolution of cluster galaxies. Recent observations have 
shown that key properties of ram-pressure stripped tails of galaxies, such as their width and structure, 
are in conflict with predictions by simulations. To increase the realism of existing simulations, we 
simulated for the first time a disk galaxy exposed face-on to a uniformly magnetized wind including 
radiative cooling and self-gravity of the gas. We find that magnetic fields have a strong effect on the 
morphology of the gas in the tail of the galaxy. While in the purely hydrodynamical case the tail is 
very clumpy, the magnetohydrodynamical (MHD) case shows very filamentary structures in the tail. 
The filaments can be strongly supported by magnetic pressure and, wherever this is the case, the 
magnetic fields vectors tend to be aligned with the filaments. Interestingly, we observe the formation 
of two dominant magnetized density tails behind the galaxy resembling the double tail observed in 
ESO 137-001. In our simulations the double tails result from the folding of the ambient magnetic 
field around the galaxy. The detectability of such structures depends on a number of factors: the 
time since the beginning of the stripping process, the characteristic length scale of the magnetic field 
fluctuations in the ICM, the orientation of the galaxy with respect to the line-of-sight, and the relative 
emissivities of the tail and the ambient ICM. Despite the fact that the effect of the magnetic field on 
the morphology of the tail is strong, magnetic draping does not suppress the rate of gas stripping, 
and the magnetic fields may in fact enhance it. Gravitational instabilities in the disk plane, combined 
with the velocity gradient between the region of the ordered field in front of the galaxy and the plane 
of the galaxy, tangle the magnetic field. In combination with the buildup of the magnetic pressure 
in front of the galaxy, this undoes the protective effect of this layer and allows the gas to leak out of 
the galaxy. The exact relationship between the amount of gas stripping in the hydrodynamical and 
MHD cases may depend on the numerical resolution and presence of an internal galactic field that is 
not included in our simulations. 
Subject headings: ram pressure stripping - galaxies: clusters: general - X-rays: galaxies: clusters 
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1. INTRODUCTION 

Galaxies populate different environments in the uni- 
verse, ranging from isolated field regions to dense galaxy 
clusters. Depending on environment, the properties of 
galaxies change: In denser regions the galaxies tend to 
contain less neutral gas, show a weaker star formation 
activity and redder colors than galaxies in underdense re- 
gions. Especially cluster spiral galaxies differ from their 
counterpa rts in the field in a number of properties (see 
[Boselli &: Gavazzi ( 2006 ) for a review) : 



• Cluster spiral galaxies tend to be redder than field 
spirals. 

• Cluster galaxies are Hl-deficient compared to their 
field counterparts. The deficiency increases to- 
wards the cluster center. Spatially resolved studies 
reveal that the HI deficiency is caused by a trun- 
cation of the gas disks. While the HI disks of field 
spirals typically extend beyond the optical disks, 
the opposite is true for Hl-deficient spirals. 



• Luminous cluster spirals have, on average, a lower 
star formation rate. The suppression of star for- 
mation goes hand in hand with HI deficiency. Spa- 
tially resolved studies reveal that the Ha disks are 
also truncated. 

• An increased 20 cm radio continuum intensity sug- 
gests an increase in magnetic field strength by fac- 
tor of 2-3 compared to field spirals. 

• Late-type galaxies follow more radial orbits and 
tend to have higher velocities than early-type 
galaxies, which suggests that they are free-falling 
into clusters. 

These observations suggest that one or more processes 
in cluster environments remove gas from galaxies or make 
them consume their gas, which leads to a subsequent 
decrease of star formation activity and change in color. 
Arguably, the most important process for this effect is 
stripping of the interstellar medium gas by the ram pres- 
sure that the galaxies experience as they move through 
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th e intracluster m edium (see also the extensive review 
by Roediger 2009 and references therein). This process 
is called rarn-pressure stripping (RPS). In order to in- 
terpret the signatures of ram-pressure stripping it is cru- 
cial to understand the fate and evolution of the stripped 
ISM. Other mechanisms that can transform spiral galax- 
ies are 'harassment' and tidal interactions. The relative 
importance of these mechanisms in transforming galaxy 
morphology, color, and gas content has been the sub- 
ject of much debate, with large sur veys helping to disen- 



tangl e the various processes (e.g. van den Bosch et al. 
2008]) . Recent deep obser vations of clusters have revealed " 



very long gas tails in HI ( Oosterloo fc van Gorkom|2005 



Koopmann et al.||2008| [Haynes et al.||2007| ). I'here have 



also been a n increasmg^umber ot observations of gas 
tails in Ha jKenney et al.'2008VGavazzi et al.'2001||Yagi 



et al. 2007 



Yoshida et al. 2008; iSun ct al. 2007) and in 



X-rays (Sun etal. 2006; Kim et~ar|2'008]|^ThrsrVikhhnin 
[2005^ Wang et al._2004J . Whether RPS can cause these 
tails and how they survive in the ICM is unknown. In 
fact, three of the long gas tails mentioned above have 
been observed i n multiple wavelengths. NGC 43 88 has a 
120 kpc HI tail (Oosterloo fc van Gorko m|2005 ), and an 
extended emission line region consisting ot many faint 
gas filaments emitting in Ha close to the galactic disc 
(Yoshida et al. 2002, 2004). Another galaxy with mul- 
tiwavelength emission is ESO 137-001 (Sun et al. 2006, 
2007), with a tail detected in both Ha and X-ray. NGC 
4438 also may have a very long tail observed in Ha (Ken- 
ney et al. 2008), while much smaller HI and CO tails 
have been observed (see VoUmer et al. 2009 and refer- 
ences therein). 

The first analyt i cal est imate of RPS appears in a paper 
by Gunn & Gott ( 1972 1 who compare the ram pressure 
with the galactic gravitational restoring force pe r unit 



area. Later hydrodynamical simulations of RPS (Abadi 
etal.|1 999, Quihs ct al. 2000, Schulz & Struck 2001 1 |Mar- 
colini ct al. 2003, .Acrcman ct al., 2003, .Rocdigcr et al. 



2006 , Rocdig cr fc ljriiggen| 120067 Roediger fc Briiggi 
2007|) confirm that this analytical estimate is fairly a 
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curate at least for galax ies t hat move prima r ily fa ce-on 



Kronberger et al. ( 2008 ) and Kapferer et al. ( 2009 1 have 
studied the impact ot RPS on star formation m spiral 



galaxi es using the SPH code GADGET. Tonnesen fc Bryan] 
( 2009 ) included radiative cooling in the adaptive hydro- 
dynamical grid code ENZD in order to investigate the im- 
pact of a clumpy, multi-phase interstellar medium. They 
modeled an inhomogeneous ISM including self-gravity 
and cooling with a resolution of 40 pc and found that 
low surface density gas is stripped easily from any radius, 
although the overall mass loss does not differ much from 
models with a smooth disc. They also find that including 
cooling results in a very different morphology for the gas 
in the tail, with a much wider range of temperatures and 
densities. The tail is significantly narrower in runs with 
radiative cooling, in agreement with observed wakes. In 
addition, they expect detectable HI and Ha signatures, 
but no observable X-ray emission for the parameters that 
they considered. 

Evidence for magnetic fields in galaxies and clusters 
comes mainly from (polarized) radio continuum emission 
and Faraday rotation measurements. Spiral galaxies are 
known to h ave a regul ar and a tangled magnetic field 
component ( Beck|[2005 and references therein) with typ- 



ical total field strengths of around lO/iC. Magnetic fields 
in the ICM are inferred from the synchrotron emission 
of cluster-wide diffuse sources in combination with in- 
verse Compton hard X-ray emission. The magnetic field 
in clusters can also be determined via the rotation mea- 
sure (RM) if the field structure and the electron densities 
are known. Magnetic fields can be studied with back- 
ground and embedded radio sources as well as with RM- 
synthesis. The ICM magnetic fields have strengths of the 
order of 1 iiG and coherence scales of the order of 10 kpc 
(e.g., |Govon i fc Feretti'2004'; 'Bon afede et al.|2010[ ). More 
recently, Kuchar & Enljlin (2011) measured the power 
spectrum in the Hydra cluster using a novel Bayesian 
method. They concluded that the magnetic field may be 
r^ 10 times stronger and detected no turnover in the mag- 
netic fluctuation power spectrum on scales from 0.3 kpc 
to 8 kpc indicating that the field correlation length is 
larger than the outer limit of that range. 

One may expect that ambient magnetic fields get 
draped around the galaxies as they move through the 
ICM. This effect of magneti c draping has been stud - 



ied by Lyutikov 



(12006) and iRuszkowski et all (12007) 



]] lp008) for the cases ot radio lobes in 
It 1 
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galaxy clusters. It "has also been pointed out that spi 
ral galaxies that move through the ICM sweep up am- 
bient magnetic field to form a dynamically im portant 
sheath arouiid them . A s pointed out, e.g., by |Pfrom-| 



mer fc Dursi] ( |2010[ ) and |Dursi fc Pfrommer| ( |2008p , the 



magnetic field strength in the boundary layer is set by 
a competition between sweeping up and slipping of field 
lines, resulting in a magnetic energy density ahead of 
the galaxy that is comparable to the ram pressure experi- 
enced by the galaxy. For typical conditions in the ICM of 
nicM ~ 10~* cm~"^ and galaxy velocities ?;gai ^ 1000 km 
s^^, this leads to a maximum field strength in the bow 

of the galaxy of J5 = \/8npv'^^^ ~ 6.5 /xG. Draping would 

modify the process of stripping, and the magnetic fields 
may prevent Kelvin-Helmholtz instabilities. The effects 
of magnetic fields on the effects of ram-pressure strip- 
ping have hardly been studied. The o nly re lated work 
is by Otmianowska-Mazur fc Vollmerl ( 2003 ) , who used 
ZEUS3U in combination with results ot sticky particle 
simulations, to study mainly the polarization properties 
of disk galaxies. Their simulation, however, neglected 
the dynamical coupling of the magnetic fields to the gas. 

2. METHODS 

2.1. Numerical techniques 

The simulations were performed with a proprietary 
version of the FLASH v3.3 code. FLASH code is a multi- 
physics parallel adaptive mesh refinement MHD code. 
Magnetic fields were comput ed using the u nsplit stag- 
gered mesh (USM) method ( |Lee fc D canc"2009; |Lee| 



20121). The USM module is based on a finite- volume 
high-order Godunov scheme (3rd order PPM) with con- 
strained transport (CT). Major strength of the USM al- 
gorithm is that it preserves vanishing divergence of the 
magnetic field to very high accuracy. We use a hybrid 
self-gravity solver. It solves for the potential on the high- 
est uniformly refined level using a parallel FFT method. 
The potential in regions refined beyond that level of re- 
finement is computed using an iterative V-cycles method. 
This hybrid approach noticeably improves code perfor- 
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mance. Radiative cooling was co mputed using Suther 
land & Dopita cooling functions (Sutherland & Dopita 
19931. We assumed solar metallicity and the cooling 
function was truncated at temperatures below 10'' K. 
Wc introduced a temperature floor at 5 x lO'^ K, and a 
density floor at 10^^^ g cm^'^. Moreover, we used an 
ideal gas equation of state with an adiabatic index of 
7 = 5/3. The mean molecular weight was allowed to 
vary as a result of changes in ionization level and was 
computed self-consistently from the Sutherland & Do- 
pita tables. Whenever the cooling timescale becomes 
shorter than the Courant timestep, we sub-cycled the 
cooling sink term assuming constant density. Although 
this method is approximate, it yields robust results and 
significantly accelerates the computations. 

2.2. Initial conditions 

In order to set up the initial distribution of gas density 
and temp erature, we follow a m ethod similar to that de- 
scribed in Melioli et al. (2008). The distribution of gas 
density is made up ot disk arid halo components. The 
halo was assumed to be isothermal with a temperature 

m n . 1 r^fiT^ . 11.-; • . 1 



Th = 8 X lO^K and a density given by 

Ph{r, z) = pm exp [-(0(r, z) - 0(0, 0))/c^ 



(1) 



where Cg is the isothermal sound speed, (/'(r, z) is the total 
gravitational potential, and pho = 2.17 x lO^^^g cm""^. 
This gas density distribution guarantees that the halo is 
in hydrostatic equilibrium. 

The gravitational potential of the galaxy consists of 
four components: dark matter halo, stellar bulge, and 
flattened stellar disk. The parameters of first three com- 
ponents are summarized in Table 1. These components 
are either analytical or are precomputed and read in at 
the code initialization stage. While the gravitational po- 
tential contribution from the gas was neglected at the ini- 
tialization stage, the gas self-gravity was included at run- 
time. In the absence of radiative cooling the gas configu- 
ration quickly settled to an equilibrium state. Apart from 
the runs corresponding to the parameter values listed in 
Table 1 and 2, we also performed two runs (one MHD and 
one hydro) with the normalization of the stellar density 
profile pd and characteristic size of the stellar compo- 
nent Vd reduced by two. These correspond to stronger 
ram pressure stripping case described below. The dark 
matter potential follows an NFW profile 
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The stellar bulge potential contribution is described by 
the Hernquist profile (Hernquist 1990) 



Mr) 



GR'h 
r + n 



(3) 



We computed the gravitational potential contribution 
corresponding to the flattened stellar disk by solving the 
Poisson equation. The density distribution in this com- 
ponent is given by the flattened King profile 
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Stellar mass distribution was truncated for (r, z) such 

that ^/{r/rdY + {z/zdY > n/rd = Zt/zd- 

The gas density and pressure of the isothermal halo 
were computed in the entire computational domain. The 
gas density in the disk was given by a sum of four compo- 
nents each of which was assumed to follow the functional 



form suggested by Wolfire et al. ( 2003 1 



p{r, z) 



2z., 



exp 



rjH 
r 



r 



(5) 



The values of rpi, r^, and Zi for all gas components (Hin, 
Hout, H2, and HI I) are given in Table 2. The transition 
between the inner Hjn and outer Hout neutral hydrogen 
density distributions was smoothed out according to p = 

g{x)pm + [1 - 9{x)]pont, where 



9{x) 



l.bx^ + 0.75x-^ 



0.25(2 




if x<0 
if < a; < 1 
if 1 < X < 2 
if a; > 2 



where x — {r — 13kpc)/lkpc. The density distribution 
given by Equation 5 was truncated using a tapering func- 
tion 



%) = 




cos(7r(r - ri)/(r2 - ri))] /2 



if r < ri 

if ri < r < r2 

if r > r2 



where ri = 20 kpc and r2 — 25 kpc. 

The boundary conditions for the density and temper- 
ature of the disk component were set to match those 
of the halo gas at Zm = 40 kpc. In order to compute 
the gas temperature in the disk, we integrated the equa- 
tion of hydrostatic equilibrium along the ^-direction from 
z — Zjnto z = kpc. This approach did not require us to 
specify the centrifugal force acting on the gas. In order 
to compute the final distribution of densities and tem- 
peratures in the entire computational domain, we either 
selected the precomputed halo gas densities and tem- 
peratures or chose the disk gas properties depending on 
which component dominated the gas pressure. The re- 
gions where the disk gas pressure dominated correspond 
to the passive scalar values equal to unity, while in the 
regions dominated by the halo pressure the passive scalar 
was set to zero. Once the gas density and temperature 
distributions have been computed, we calculated the gas 
velocity required to keep the gas in rotational hydrostatic 
equilibrium from 



rdP 
p dr ' 



(6) 



where v^ — rd(t)/dr. While the disk component of the 
gas follows this velocity curve, the halo component was 
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TABLE 1 

Potential parameters 



Pdm 

[10-24 g cm-3] 


Pd 

[10-24 g cm-3] 


A/fc 
[IQIOMq] 


l-dm 

[kpc] 


rb 
[kpc] 


rt 
[kpc] 


rd 
[kpc] 


[kpc] 


Zd 

[kpc] 


2.9 X 10-1 


1.2 X 102 


3.5 


3.08 X 10^ 


0.8 


6.44 X 10^ 


4.83 


1.93 


1.16 X 10-1 



TABLE 2 
Density distribution 



Component 


Ed[MQ/pc2] 


roi [kpc] 


r,[kpc] 


^.[Pc] 


Hill 


2.38 X IQi 


1.0 


l.Ox^ 


1.78 X 102 


Hout 


1.71 X 10^ 


1.0 X lOi 


4.0 


3.24 X 102 


H2 


1.72 X 102 


3.3 


2.89 


6.34 X 10^ 


HII 


4.17 


0.0 


3.0 X 10^ 


8.80 X 102 



assumed to be initially static. IVIodels corresponding to 
the values listed in table 1 and 2 have disk gas mass of 
~ 3 X 10^° 1VI0, stellar mass ^ xlO^^ IVI©, and halo mass 
~ xlOi2 Mq. 

The computational domain was 120 kpcxl20 kpcx240 
kpc. The disk of the galaxy was parallel to the x — y 
plane and its center was located at z = 50 kpc. We 
used static mesh refinement. The domain was fully re- 
fined for z < 60 kpc and had one level of refinement less 
throughout the rest of the volume. This ensured that 
the disk was always maximally refined. It also allowed 
us to take advantage of the fast gravity solver in the most 
of the volume providing noticeable improvement in the 
code performance. The maximum refinement level was 
6 (low resolution) and 7 (high resolution). These refine- 
ment levels correspond to maximum spatial resolutions 
of 469 pc and 234 pc, respectively, and the first zone cen- 
ter is half the resolution element size away from the disk 
midplane. We compared the disk gas mass (defined as 
having passive scalar equal to unity) in the initial state 
for our two resolutions and found that the masses differed 
by less than 3%. While the refinement blocks have the 
same aspect ratio as the whole computational domain, 
each zone is cube-shaped. We note that in the adaptive- 
mesh refinement tests the entire domain got refined in 
just one shock crossing time. In that mode, the use of 
the fast hybrid gravity solver was not warranted as the 
gravity solver was accelerated only on a small number 
of blocks that had a low refinement level. Therefore, we 
decided to use highly refined static mesh in conjunction 
with the hybrid gravity solver. 

The equilibrium gas density, temperature, and velocity 
were precomputcd on a grid with 25 pc resolution in the 
z-direction and 50 pc in the radial direction. In order 
to set up the initial state, we used a gravitational po- 
tential precomputcd on a logarithmic grid over a range 
from 15.1 pc to 515.8 kpc in both r and z using 400 
logarithmically spaced points in each direction. These 
initial conditions were mapped onto the simulation grid 
at the code initialization stage. The tabulated values of 
the gravitational acceleration due to the static potential 
were also computed in advance at very high resolution. 
These tables were then used to interpolate the gravita- 
tional acceleration at runtime. In addition to the acceler- 
ation due to the static potential, the gas was also allowed 
to self-gravitate. We verified that the initial conditions 
were globally stable, and the system quickly settled to 



an equilibrium state over a few orbital times in the ab- 
sence of radiative cooling. While we experimented with 
disks that have an initial Toomre Q parameters greater 
than unity, the parameters in Table 1 and 2 correspond 
to gravitationally unstable disks that fragment and in- 
crease their porosity over time. 

2.3. Boundary conditions 

We used isolated boundary conditions for the 
hybrid multigrid-|-PFFT gravity solver (using 
Dirichlet boundary for the multigrid). That is, 
we set grav_boundary_type to "isolated" and 
mg_boundary_type for all boundaries to "dirichlet" 
in the FLASH parameter file. Outflow boundary con- 
ditions for the gas are used on all but one boundary. 
On the inflow boundary we assumed that the gas had 
the same temperature as the halo gas, T = 8 x 10^ 
K, a density of Pwind = 1-67 x 10~^^ g cm-'^, and a 
maximum inflow speed of Wmax = lO^cm s~^. This 
yields a IVIach number of inflow of M = v/cg ~ 2.3, a 
plasma (3 = p/{B'^/87r) w 3.6 and an Alfvenic IVIach 
number of Ma = My/(ffj2 « 4.0. Note that the ratio 
of the wind ram pressure to the magnetic pressure 
^ram/^i3 = 2Af^ « 31.5 » 1. The inflow speed 
Vhi{t) = finit)vnia.x was increased gradually from to 
lO'^km s~^ as a function of time t according to the 
following profile 



/in(t) = 1 



l-1.5a;^-^0.75x^ 
0.25(2 -a;)3 




if a; < 1 

if 1< a; < 2 

if a;> 2 



where x = t/At and At = 58.6 JVIyr. The direction of 
the inflow was along the z-direction and parallel to the 
long side of the computational box. The magnetic field 
injected at the boundary was in the x — y plane, such that 
Bx/V^ = By/\/4n = 10~^/a/2 Gauss. This choice was 
motivated primarily by simplicity. The actual topology 
of the magnetic field in real galaxies is either not known, 
or its description would require us to introduce addi- 
tional free parameters to the model. The morphology 
of the internal steady-state field is likely to be governed 
by a galactic dynamo, the study of which is beyond the 
scope of our current work. The simulated steady-state 
internal magnetic fields may also depend on other fac- 
tors (e.g., the presence of winds) that affect the relative 
strength of the global toroidal and poloidal galactic field 
components in a differentially rotating galaxy. Thus, we 
decided to start with the simplest setup in which the ef- 
fect of the external field s can be simulated. A si milar 
approach was adopted by Pfrommer & Dursi ( 2010| ) who 
considered only external fields but approxmiated galax- 
ies as nonrotating high density clouds. We extend their 
work by including galactic rotation, hot galactic halo, 
radiative cooling, and self-gravity. We defer the investi- 
gation of the impact of internal galactic fields on the ram 
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to unity may underestimate the volume occupied by the 
galactic gas, which in turn would lead to the overesti- 
mate of the amount of gas stripping. We also note that 
Tonnesen & Bryan (2009) use the same threshold for the 
passive scalar to identify the ram pressure stripping tails 
in their simulations, and also conclude that the results 
remain unchanged if a threshold of 0.1 is used instead. 

The rebound in Mgas around 0.25 Gyr is due to the 
gas fallback after the initial episode of ram pressure 
stripping (see also Roediger et al. 20061. In general, 
the hydro case exhibits a relatively niodest amount 

_2g - of gas stripping, which is not unexpected. Following 

Gunn & Gott (1972), we can estimate the gravita- 
tional restoring force per unit surface acting o the 

-26.'^ S^^ -1 grav = 27rGzjgas^starsj wherc ijgas and X/stars are 

the surface density of the galactic gas and the surface 
density of stars, respectively. We can then compare 
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Fig. 1. — Gas density projected along the line of sight. Left: 
MHD case. Right: non-magnetic case. Both cases are for t = 0.8 
Gyr. The difference between tail morphologies in these two cases 
is striking. 

pressure stripping process to a future publication. 

The simulations were performed on the Pleiades ma- 
chine at NASA Ames on 1024 processors. Higher resolu- 
tion hydro run required approximately ^ 10^ CPU hours 
to advance to 0.9 Gyr, while its MHD counterpart took 
roughly ~ 3 x 10^ CPU hours to reach 0.9 Gyr. 

3. RESULTS 

3.L Density tail morphology and gas stripping rates in 
MHD and hydro runs 

Figure IT] shows the gas density projected along the line 
of sight (Indirection) 0.8 Gyr after the start of the wind. 
The left panel corresponds to the MHD case and the 
right panel to the non-magnetic case (all images in this 
paper correspond to high resolution runs). The difference 
between tail morphologies in these two cases is striking. 
While the purely hydrodynamical case the tail is very 
clumpy, the MHD case shows very filamentary structures. 
As a result of this large-scale inhomogeneity in the MHD 
case, the surface density shows more variation with the 
orientation of projection in the sense that the number of 
filaments seen in projection depends on the line of sight. 

Figure |2] quantifies the amount of ram pressure 
stripping of the galactic gas. Shown as a function of 
time is the disk gas mass normalized to its initial value. 
The mass in the disk is computed within a rectangular 
solid of width and depth equal to 50 kpc, and height 
equal to 6 kpc, that is centered on and parallel to the 
disk, and where the the galactic passive scalar values are 
greater than 0.25. The disk gas mass is then obtained 
by integrating the gas density multiplied by the passive 
scalar values in this volume. Thick lines correspond 
to the higher resolution cases. Solid and dashed lines 
are for the MHD and hydro cases, respectively. In the 
hydrodynamical case, this prescription for calculating 
the retained mass leads to results that are entirely 
independent of the passive fluid threshold as long as it is 
less than ^0.5. The MHD case is even less sensitive to 
this threshold. Wc point out that threshold values close 
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instance, near r ~ 13 kpc in the plane of the disk where 
the transition between two neutral hydrogen density 
regimes Hin, Hout occurs, Pgrav exceeds Pram by about an 
order of magnitude. Therefore, the amount of stripping 
is expected to be moderate. For comparison, we also 
computed two cases where stripping is stronger (as 
mentioned in Section 2.2). These runs were computed 
at the same resolution as our low resolution runs. 
Specifically, we performed one hydro and one MHD 
run for a weaker stellar gravitational field, comparable 
magnitude of ram pressure, and weaker magnetic field. 
These runs had 3 times smaller wind magnetic field and 
density, and v3 faster wind. The normalization of the 
stellar mass profile pd and the characteristic radius of 
the disk r^ were reduced by 2. These runs resulted in 
stronger gas removal (see right panel in Figure ^ . 

The most interesting point from both panels in 
Figure [2] is that the MHD case shows larger amounts of 
gas stripping than the purely hydrodynamical case. This 
suggests that the presence of magnetic fields accelerates 
the removal of gas, contrary to the expectation that 
magnetic draping reduces the amount of ram pressure 
stripping. 



3.2. Formation of the ordered field ahead of the galaxy 
and filamentary field in the tail 

In Figure p^ we present cross sections {y — 0) through 
the computational domain showing the absolute value 
of B, (left panel) and B^ at t ^ 0.45 Gyr. This fig- 
ure clearly shows the detached shock that formed after 
about 0.2 Gyr. The jump in the z-component of the gas 
velocity across the shock measured along the galaxy's ro- 
tation axis at i = 0.45 Gyr is ^ 2.3. A draping layer is 
established closer to the disk plane than the shock stand- 
off distance. The thickness of the draping layer may be 
defined as the width of the region where the sum of the 
gas pressure and magnetic pressures become comparable 
to the ram pressure (Pram ^ 1-7 x 10""'^"'^). The vertical 
magnetic field reveals the presence of elongated filaments 
that extend to large distances from the disk of the galaxy. 

In Figure |4] we show the profiles of the magnetic pres- 
sure Pb across the galactic disk plane a,t t — 0.45 Gyr. 
The magnetic pressure is computed within cylinders of 
radii 25 kpc (dotted), 15 kpc (solid), and 5 kpc (dashed). 
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Fig. 2. — Left: time dependence of the galactic disk gas mass normalized to its initial value. Solid and dashed lines are for the MHD 
and hydro cases, respectively. Thick lines correspond to the higher resolution cases. Right: same as left panel but for low resolution strong 
stripping case (see text for details). 
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Fig. 3. — Cross sections (y = 0) through the computational 
domain showing the absolute value of the vertical magnetic field 
Sz/v47r (left) and the horizontal field Bx/V^^ in /^G. Both 
panels are for t = 0.45 Gyr. 

Thin horizontal dashed hne corresponds to O.SpwindW^ind- 
The thickness of the draping layer corresponds to the 
point where Pb = O.SpwindwJind- This definition is comes 
from the assumption of the equipartition between the 
magnetic and thermal pressures in the layer. This as- 
sumption is better when radiative cooling is slow and we 
verified that it approximately holds in the draping layer 
in front of the galactic disk plane. 

The thickness of the layer decreases with the distance 
from the rotation axis of the galaxy. For a spherical ob- 
stacle of radius r, the theoretical thickness of the drap- 
ing layer along the line connect ing the observer and the 



stagnatioii poiiit is I ^ R/M'j ( Lyutikov 



Pfrommer 2008 



2006 



Dursi fc 

Pfrommer & Dursi||201U|). However, the 
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geometry of the problem is ditterent m our case. Our thin 
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Fig. 4. — Profiles of the magnetic pressure Pg across the 
galactic disk plane. Front side of the galaxy corresponds to the 
positive values of the distance from the disk midplane. P^ is 
computed within cylinders 25 kpc (dotted), 15 kpc (solid), and 5 
kpc (dashed) in radius. Thin horizontal dashed line corresponds 
to 0.5p„ind''^inti ~ 0.5Pram. The thickness of the draping layer 
corresponds to the point where Pb = 0.5Pram. 

disk is oriented face-on with respect to the wind. A finite 
disk may be better approximated by a part of a sphere of 
very large curvature radius R > r, and so the thickness of 
the draping layer is much larger in the vicinity of a disk 
than close to a spherical obstacle of radius r. For a finite 
disk, the central thickness of the ma gnetic layer can be 
approximated following .Erkaev et al. (1995J . From mag- 
netic freezing, and assuming one-dimensional flow near 
the flat disk, we get i^wind'^wind = B{zd)v{zd), where Zd is 
the distance from the galactic plane, and v is the gas ve- 
locity normal to the disk plane. It can be shown that the 
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Fig. 5. — Line-of-sight projections of the ratios of the magnetic and thermal pressures. Result for passive scalar fraction less than 0.5 
(more than 0.5) is shown in left (middle) panel. Both cases are for t = 0.8 Gyr. Right: zoom-in on the regions immediately above the 
disk plane shown in the middle panel Magnetic field vectors are superimposed on the image and their length is proportional to the field 
strength. 



velocity linearly approaches zero as a function of the dis- 
tance from any su rface along the stagnation line ( Landau 
fc Lifshit"z~||1978 ). Specifically, v = a/xuwind, where a is a 
dimensional constant around unity that weakly depends 
on the surface geometry (e.g., a = 3 for a sphere, a = 2 
for a cylinder; assumed equal to unity below), and /i is 
the dimensionless distance (ratio of the distance to the 
characteristic size of the obstacle). Thus, the magnetic 
field is B{zg) — Bwind//^- Comparing magnetic pressure 
to half the ram pressure, and assuming equipartition and 
neglecting cooling, we get 



R2 1 

^wind _ ^ 



Sna^fi 



2 ~ r)Pwind''^wind 



(7) 



The dimensionless thickness of the layer is thus fi = 
MJ , which is much larger than the estimate for a 
sphere. For a characteristic size R ^ 25 kpc of the disk, 
the central thickness is Z ~ R/Ma ~ 6 kpc. Close to the 
edge of the disk the thickness of the layer is expected to 
be smaller than this rough estimate due to the smaller 
effective curvature radius of the surface and the slippage 
of the magnetic field lines past the disk. These trends are 
supported by the results shown in Figure |4] that clearly 
shows that the typical thickness is around ~ 5 kpc (mea- 
sured within 15 kpc radius cylinder; solid line) which 
is comparable to our estimate. As expected, the layer 
thickness decreases away from the rotation axis, and is 
easily resolved in our simulations. The magnetic layer 
forms on a timescale iform ^ V'^^wind ^ ^/(^'wind-^A) ^ 6 
Myr. In the low-resolution runs for the strong wind case 
(see right panel in Figure [2]) the central thickness of the 
draping layer was 3 times smaller than in our fiducial 
case due to weaker magnetic field and faster wind, but 
was nevertheless still resolved. 

The velocity gradient between the disk and the ICM, 
tangles the field close to the disk midplane and below the 
draping layer. This tangling is clearly seen in both the 



Bz (left panel in Figure pi) and B^ maps and it is caused 
by shear instabilities whose growth rate is given by 



27r ^ 

WKH = -—i\v — 

A Pi, 



Pic 



isniA^icm 



Pu 



(8) 



where A is the wavelength of the unstable mode and Aw 
the velocity difference across the shear layer. Taking 
Pism ^ Picmi wc Can Write for the e-folding growth time 



shear 



^KH 



40Myr 



A 200km/s ^pism/Pi< 
1 kpc Av 50 



,(9) 



This timescale estimate is more appropriate after the 
wind hits the galaxy. At earlier times the density con- 
trast between the ISM and ICM is larger and the ve- 
locity shear between the dense disk and the lower den- 
sity co-rotating ambient gas is considerably smaller than 
the fiducial value used in the above estimate. We note 
that the shearing layer may in principle be stabilized to 
some extent by the disk gravity. However, the Richard- 
son number Ri that measures the ratio of potential to 
kinetic energy is 



Ri 



gzh 
Av2 



(10) 



where h is the disk thickness, g^ is the vertical compo- 
nent of the gravitational field in the disk, and where we 
assumed that the shearing velocity Aw is comparable to 
the orbital velocity. Since the disk is thin, Ri <^ 1 and 
gravity will not prevent the KH instability. However, the 
KHI can also be suppressed by the finite thickness of the 
shearing layer. General arguments suggest that modes 
kd > 1, where d is the thickness of the shearing layer, 
are suppressed ( Chandrasekhar|1961 Ferrari et al.|1982" 
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Fig. 6. — Volume rendering of the magnetic field strength and the gas density. The disk of the galaxy is seen edge-on in the left column 
and face-on in the right one. Top row corresponds to the magnetic field strength and the bottom one to the density. The snapshot is taken 
at 0.45 Gyr. The seven iso-surfaces in these volume rendering plots are distributed uniformly on a logarithmic scale of base 10 from -24 
to -28.5 (for the gas density) and linearly from 6 X 10~^ to 1.5 X 10"*^ GsiUss/\/4n (for the magnetic field strength). Double tail in the 
magnetic field strength and density is clearly seen in both panels. 



Churazov & Inogamov 2004). For d ~ 1 kpc, pertur- 
bations characterized by A ;^ a few kpc should be sup- 
pressed. This sets a lower limit on A in Eq. (9). We also 
point out that the growth of the KH unstable modes will 
be limited by the onset of nonlinear saturation, and the 
modes much shorter than the disk thickness are unlikely 
to significantly perturb the disk. Magnetic fields will 
suppress unstable modes parallel to the field but will not 
suppress modes perpendicular to the field lines. However, 
in our setup the disk galaxy shears against the ambient 
ICM and the magnetic field will not be able to suppress 
the onset of instabilities. 

The galactic disk can also be unstable to the Rayleigh- 
Taylor instability (RTI). The infalling galaxy is subject 
to the ICM drag force 



^d — rj^'P^wind'^'wind^i 



(11) 



where Cp is the cofScient of parasite drag, P wind = Pi 



and A is the cross-s ection of the obstacle (von Mises 
1959 Batchelor||2000 )■ For a circular disk oriented tace- 
on with respect to the wind Cp = 1.11. This drag force 
could cause the RTI. However, as po i nted o ut by sev- 
eral authors (e.g., Roediger & Hensler (2008) or Murray 
et al. (1993)), the gravitational acceleration gdisk of the 
disk can suppress the RTI if gdiskAhp > i^drag, which 
leads to the same condition as that for the absence of 



the ram pressure stripping, namely .gdisk^disk ^ ^ram- 
Thus, the gas does not become RT unstable before it 
gets stripped. Even the stripped gas may potentially be 
stabilized against RTI by the self-gravity of the gas if 



27rl], 



disk 



> 



Cpp 



pPicm'^wind 



(12) 



This implies that the stripped gas is stable against RTI 
if the surface density exceeds Scrit given by 



Cp J ram 

4^ G 



1/2 



(13) 



The value of Scrit is comparable to Sdisk in our simu- 
lations. Finally, magnetic fields may also stabilize the 
growth of the perturbations along the field lines if the 
perturbations c orrespond to wavelen gths smaller than a 
critical length ( Chandrasekhar 1961). However, pertur- 
bations in the direction perpendicular to the magnetic 
field will not be suppressed. Three-dimensional simula- 
tions demonstrate that the overall RTI growth rate will 
not be supressed and may in fact be somewhat acceler- 
ated in the MHD case ( [Stone fc Gardiner||2007 1 . 

The timescale of the gravitational instability, igrav , is 
comparable to the orbital timescale iorb — ''/^gai (orbital 
period is ^ 600 Myr at r = 20 kpc). We also point out 
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Fig. 7. — Same as Figure|6]but for t = 0.35 Gyr. 



that the formation timescale iform of the draping layer 
is shorter than all of the above timescales, so the hierar- 
chy of timescales is: tform < ishear ^ igrav All of these 
timescales are significantly shorter than the simulation 
time. The estimates presented above are consistent with 
the observation that no clear disk fragmentation occurs 
before ^ 50 to 100 Myr when the wind reaches the disk, 
and with the fact that the disk is clearly unstable at 
~ 0.35 Gyr. 

The build-up of the inhomogeneous magnetic field 
causes pressures that undo the protective effect of the 
draping layer and allow the gas to stream out of the 
galaxy. In the process, the ordered horizontal fields are 
transformed into vertical filamentary fields that fiuctu- 
ate on small scales. Small-scale vertical fields are seen 
very close to the disk plane in the left panel in Figure [3J 
In the downstream region, vertical fields are formed that 
stretch over large distances from the galactic disk. We 
speculate that the gas removal may be further enhanced 
because the magnetic fields anchored in the tail may pull 
the gas out of the galaxy. The strong amplification of 
the magnetic field in the disk plane is clearly seen in 
Figure [4] This is caused by the combination of the flux 
freezing, gas compression due to ram pressure of the su- 
personically moving galaxy, radiative cooling, and the 
differential motions. The increase in the gas density is 
due to the ram pressure of the wind and the radiative 
cooling which allows for the gas condensation in the disk 
plane. 

Interestingly, the trends with numerical resolution dif- 
fer in the hydro and MHD cases. In the purely hydro- 



dynamical case, the efficiency of stripping increases with 
resolution, while the MHD case shows the opposite be- 
havior. Why the amount of stripping becomes less with 
increasing resolution in the MHD case is not clear. Runs 
at even higher resolution are currently unfeasible. It is 
due to this limitation that we cannot draw firm conclu- 
sion on the relative magnitude of the ram pressure strip- 
ping in the hydrodynamical and MHD cases. Also, we 
should note that we are using a simplified setup in which 
the ICM is threaded by a uniform magnetic field, and we 
ignored any internal galactic magnetic fields. Such inter- 
nal fields may tend to protect the galaxy against strip- 
ping. However, even though our setup does not include a 
magnetic field internal to the galaxy, the ambient large- 
scale external field gets amplified in front of the galaxy 
over time, and yet the stripping is not inhibited. We 
conclude that further analysis for a range of wind and 
galaxy parameters, and accompanied by a careful con- 
vergence study (e.g., in smaller computational volumes), 
is needed to test our hypothesis that the magnetic field 
may accelerate the gas removal, or at least not necessar- 
ily impede it. In an upcoming foUowup study we will 
present observational diagnostics - HI, Ha, X-ray signa- 
tures - for a range of the parameters characterizing the 
wind and galaxy, including disk orientation and intrinsic 
galactic fields. 

3.3. Magnetic pressure support and field orientation in 
the tail filaments 

Figure |5] shows the line-of-sight projections of the ra- 
tios of the magnetic and thermal pressures at 0.8 Gyr. 



10 



Ruszkowski et al. 



Left panel corresponds to the contribution from the gas 
associated with the passive scalar fraction below 0.5, and 
the middle one corresponds to the passive scalar values in 
excess of 0.5 (i.e., stripped gas). This figure demonstrates 
that the filaments in the tail can be strongly supported 
by the magnetic pressure. A zoom-in on the regions im- 
mediately above the disk plane shown in the right panel. 
This enlarged region reveals that the magnetic field vec- 
tors (shown as arrows overlaid on the image; arrow length 
proportional to the field strength) tend to be aligned with 
those regions where the magnetic pressure support dom- 
inates. This hints at the interesting possibility that this 
effect could prevent anisotropic thermal conduction from 
evaporating the magnetically supported filaments. 

3.4. Formation of double tails 

In Figure [6] we show the volume rendering of the mag- 
netic field strength (top panels) and the gas density. In 
the left column, the disk of the galaxy is seen edge-on, 
and the right column corresponds to the face-on orienta- 
tion of the disk. The snapshot is taken at 0.45 Gyr and 
corresponds to the high-resolution runs. The seven iso- 
surfaces in these volume rendering plots are distributed 
uniformly on a logarithmic scale of base 10 from -24 to - 
28.5 (for the gas density) and linearly from B/^/At: — 1.5 
/zGauss to B/^/in = 6 ^Gauss (for the magnetic field 
strength). The upper left panel in this figure clearly illus- 
trates that the magnetic field is ordered on the working 
surface of the galactic disk (i.e., in front of the galactic 
disk) and that the field gets tangled immediately behind 
the disk and then stretched along the direction of the 
wind. Two magnetic tails are clearly seen behind the 
galaxy in the left column panels. Moreover, there is a 
tendency for the gas to spatially correlate with the mag- 
netic tails. While the gas distribution is in general fila- 
mentary, it is clear that two dominant density tails are 
present, and that they are denser than the ambient ma- 
terial. The density contrast between the ambient ICM 
and the tails measured near the end of the tails is up 
to 40%. We note that the contrast could certainly be 
considerably higher in the simulations characterized by 
lower ICM wind density but with the magnitude of Pram 
still comparable to the one adopted here (and assuming 
comparable wind gas pressure and magnetic field). This 
suggests that the stripping process could generate an ob- 
servable double tail depending on the orientation of the 
galaxy's direction of motion with respect to the linc-of- 
sight. Such a double tail behind a galaxy has indeed been 
observed in ESQ 137-001 (Sun et al. 2006, 2007). In or- 
der for the tails to be detectable in X-ray or Ha data, one 
needs to consider a range of wind and galaxy parameters 
to assess if the tails can be bright enough to be detected, 
and that there is sufficient contrast between the tails and 
the ambient ICM emission. This task requires additional 
(and very CPU-intensive) runs and we defer this study 
to a follow-up paper. Here we merely point out that the 
dynamical interaction of the wind and the galaxy may 
naturally lead to the formation of bifurcating tails if the 
external magnetic field is uniform on scales sufficiently 
large compared to the disk size. We also note that for 
galaxies oriented edge-on with respect to the wind, there 
exists an asymmetry between the side of the galaxy mov- 
ing upstream and downstream relative to the wind. This 



may further complicate the tail properties. However, the 
rotational speed of the galaxy is substantially smaller 
than the infall speed of the galaxy in the ICM. 

In order to understand the morphology of the tail, we 
compared snapshots taken at 0.35 Gyr and 0.45 Gyr. 
Snapshot taken at 0.35 Gyr is shown in Figure [7j Given 
the wind speed of 1000 km/s, the time difference of 0.1 
Gyr between these snapshots translates to 43% of the 
longer box length Lz {L^ = 240 kpc). At the same time, 
the typical rotational velocity of the galaxy is around 
^ 200 km/s, so in 0.1 Gyr the galaxy rotates by ~ 16% 
of the circumference at r = 20 kpc. A comparison of 
these two snapshots reveals that the endpoints of the 
tails remain on the diagonal line in the x — y plane at 
z — Lz- This diagonal line coincides with the orientation 
of the external magnetic field in the wind. This shows 
that the two magnetic tails are due to the folding of the 
external field around the galaxy. Thus, this could be a 
generic feature whenever the coherence length of the ex- 
ternal field is larger than the characteristic radius of the 
galaxy. 

We emphasize that the tails are not due to "limb 
brightening" of a hollow cone of the gas density or mag- 
netic field distribution seen in projection onto the plane 
of the sky. Face-on projections clearly show that the 
tails lie along the diagonal plane diagonal x — y plane 
(see right column panels in Figure [6|. However, at later 
times the two distinct tails get washed out or they are 
not as cleanly separated. The expectation is that this 
happens in the time it takes for the disk edge to ro- 
tate by TT radians, i.e., in about r^ 0.3 Gyr since a clear 
double tail has been established at ~ 45 Gyr. In fact, 
snapshots taken around ~ 0.8 Gyr do not reveal the dou- 
ble tailjwhich is consistent with the above estimate (see 
Figure^. However, the tail section of the flow is nev- 
ertheless still somewhat flattened along the diagonal line 
in the x — y plane. 

3.5. Phase space distribution for the gas and the 
magnetic field 

The left panel in Figure [9] shows a B — p phase space 
for the gas downstream from the galaxy at 0.8 Gyr. The 
color bar shows the amount of gas in solar masses in 
each B ~ p bin. The right panel shows an analogous 
quantity for the gas associated with the disk fluid frac- 
tion greater than 0.5 (i.e., stripped gas). These phase 
space diagrams demonstrate that: (1) there are no high 
density— low magnetic held strength regions, and (2) the 
gas in the tails tends to be slightly more magnetized than 
the ambient ICM. Comparison of both panels reveals that 
the right maximum in the B — p plane corresponds to the 
galactic gas. That peak in the distribution corresponds 
to slightly but systematically higher magnetic fields than 
the left peak that is attributed to the ambient gas. Sim- 
ilar two-dimensional histograms are shown in Figure |10[ 
Both panels in this figure show T—p. The left panel is for 
the MHD case and the right one for the hydro case. The 
phase diagrams show that the gas is clearly multiphase. 
Significant amounts of gas near T ^ lO^K are due to the 
isothermal halo and the wind. Evidence for adiabatic gas 
compression is clearly seen. Note that the MHD case has 
more pockets of low density-high temperature gas. 

4. CONCLUSIONS 
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Fig. 8. — Volume rendering of the magnetic field strength (top 
panel) and the gas density at t = 0.8 Gyr. The disk of the galaxy 
is seen edge-on. Magnetic field and density contours are at the 
same levels as in Figure [6] and Figure [7] 

We performed hydrodynamical and magnetohydrody- 
namical simulations of ram pressure stripping in a disk 
galaxy including radiative cooling and self-gravity of the 
ISM. Our main conclusions can be summarized as fol- 
lows: 

1. The presence of ambient magnetic fields has a 
strong impact on the morphology of the tail. While 
the MHD case shows long (> 100 kpc) filamentary 
structures in the tail, in the purely hydrodynami- 
cal case the tail is very clumpy on scales of a few 
kpc. 

2. In the filaments of the ram-pressure stripped tail 
the magnetic pressure support dominates and is 
aligned with the filaments. This may prevent ther- 
mal conduction between filaments and the ambi- 
ent, hot ICM and could be relevant to the question 
whether stars can form in ram-pressure stripped 
gas. The filamentary structure helps explain obser- 
vations of ram-pressure stripped tails. Our setup 
is not suitable to reproduce specific observations. 
However, our finding that the galactic wake con- 
sists of long filaments that are magnetically sup- 
ported helps to explain observations such as those 
by Yoshida et al. (2008) who find narrow, blue fil- 
aments with Ha emission that extend up to 80 kpc 
from a galaxy in the Coma cluster. 



3. The ram pressure stripping process may lead to 
the formation of double, magnetized density tails 
behind the galaxy. In our simulations the bifur- 
cated tails result from the folding of the ambient 
magnetic field around the galaxy. The detectability 
of such structures depends on a number of factors: 
the time since the beginning of the stripping pro- 
cess, the characteristic lengthscale of the magnetic 
field fiuctuations in the ICM, the orientation of the 
galaxy with respect to the line-of-sight, and the rel- 
ative emissivities of the tail and the ambient ICM. 



4. The presence of uniform magnetic fields in the ICM 
leads to larger amounts of gas stripping than in 
the absence of magnetic fields (for the numerical 
resolutions considered here). A disk galaxy shears 
against the ambient ICM and the magnetic field 
will not be able to suppress the onset of insta- 
bilities. Even if the headwind experienced by the 
galaxy is attenuated by the magnetic draping layer 
(whose interface is stable) , the pressure on the disk 
suffices to remove the gas from the galaxy. In 
the process, the ordered horizontal fields are trans- 
formed into vertical filamentary fields that fiuctu- 
ate on small scales. We speculate that the gas re- 
moval may be further enhanced because the mag- 
netic fields anchored in the tail may pull the gas 
out of the galaxy and the magnetic forces acting 
on the edge of the galaxtic disk excert additional 
forces that may help to remove the gas. Hence, 
magnetic draping may not reduce the amount of 
ram pressure stripping. This situation may be dif- 
ferent in giant elliptical galaxies where the rotation 
does not play a significant role. Higher resolution 
simulations are necessary to test this hypothesis. 
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